AN ALGORITHM FOR COMPUTING BOUNDARY EXTENSIONS 

OF CONFORMAL MAPS 
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Abstract. We show that from a sufficiently good approximation to a con- 
formal map of the unit disk onto a bounded domain D, a sufficiently good 
approximation to the boundary of D, and sufficient local connectivity infor- 
mation for dD, it is possible to compute arbitrarily good approximations to 
the boundary extension of (p. 



1. Introduction 

Let C denote the complex plane, and let D denote the disk whose center is the 
origin and whose radius is 1. A domain is a subset of C that is open and connected. 
Let N denote the set of natural numbers {0, 1,2,.. .}. 

A compact set X C C is uniformly locally arcwise connected if there is a function 
5 : N — ?> N such that whenever p, q are distinct points in X and |p — g| < 2^^('^\ 
X contains an arc A from p to g whose diameter is smaller than 2~^. When X is 
compact and connected, this is equivalent to local connectivity. See e.g. Chapter 3 
of [5]. We call the function g a ULAC function for X. 

Suppose D is a bounded and simply connected domain. By the Riemann Map- 
ping Theorem, there is a conformal map of ID onto D, (p. If the boundary of D 
is locally connected, then ip extends to a continuous map of D onto D. We call 
this map the boundary extension of (f) and denote it by (j) as well. It follows from 
the converse of the Hahn-Mazurkiewicz Theorem that (f) has a boundary extension 
precisely when the boundary of D is locally connected. If the boundary of D is a 
Jordan curve, then the boundary extension of is a homeomorphism. This is the 
famous Caratheodory Theorem. These results are proven in jT3. . 

Boundary extensions of conformal maps play a key role in the solution of Dirichlet 
problems for the domain D. It is thus natural to investigate possibilities for their 
computation. To this end, let us discuss how we shall approximate points, compact 
sets, and continuous functions from D into C. 

A rational rectangle is a rectangle whose vertexes are rational points. Let us agree 
to approximate points in C by rational rectangles. That is, a rational rectangle R 
approximates a point p G C if p G R. 

A finite cover {Ui, . . . ,Uk} of a compact set X by open sets is tight if each Uj 
contains a point of X. Let us agree to approximate compact sets by tight covers 
of rational rectangles. Thus, sufficiently accurate approximations to a compact set 
allow one to plot it on a computer screen with arbitrarily high resolution. 
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Whenever < ri < 7-2 < 1 and 9i < 62 are such that 62 — 61 < 27r, let 
Cri,r2,ei,e2 = {re'^ ■ ri < r < ra A Oi < 9 < 6*2}. 

We call such a set a Carleson rectangle. If ri, r2, 0i, 02 arc rational, then we call 
it a rational Carleson rectangle. 

A Carleson rectangle does not contain 0. So, let us use the term rational disk 
to refer to an open disk whose center is a rational point and whose radius is a 
rational number. We could approximate points by rational disks. However, rational 
rectangles have certain useful properties such that the intersection of two rational 
rectangles is a rational rectangle. 

A continuous function / : D — > C can thus be approximated by a finite set of pairs 
{{Ui,Vi), . . . , ([/„, Vn)} such that each Uj is either a rational Carleson rectangle or 
a rational disk whose center is the origin, Uj C D, and such that f[Uj] C Vj for 
each j. Rational rectangles could be used instead of Carleson rectangles; the choice 
is for the sake of convenience later. 

A key consequence of these approximation schemes is that from sufficiently good 
approximations to p and / one may compute arbitrarily good approximations to 
f{p) in a fairly direct manner. A more theoretical case for these approximation 
schemes is given in [3] and |14) . 

Finally, let us agree to approximate a function 5 : N — >■ N by finite initial 
segments of g. That is, a restriction of g to a finite initial segment of N. 

We can state our main result as follows. There is an algorithm that, from ap- 
proximations to (j), dD, a point p G D, and an approximation to a ULAC function 
g for dD, computes an approximation to 4'{p). Furthermore, the algorithm con- 
verges in that it yields arbitrarily good approximations to 4>{p) when provided with 
sufficiently good approximations to (j), dD, p, and g. 

This appears to be the first general algorithm for computing boundary extensions 
of conformal maps. Hence, the results herein evidence that it is possible to create 
software for computing boundary extensions of conformal maps. We note that 
no amount of differentiability of the domain boundaries is assumed; all that is 
required is sufficient information about the local connectivity of the boundary. In 
[7], P. Hertling demonstrated the computability of the Riemann mapping operator 
for simply connected domains. The proof is based on Koebe's construction and 
the work of Bishop and Bridges [9], [2]. Algorithms such as those in [10] and [1] 
compute approximations to conformal maps of the unit disk onto Jordan domains 
from approximations to the domain boundary. In , it is shown that the boundary 
extensions of such maps can not be computed from approximations to the boundary 
alone, and moreover that a ULAC function for the boundary represents both the 
sufficient and necessary additional information required. In particular, there is a 
Jordan domain D whose boundary is computable but such that no conformal map of 
D onto D has a computable boundary extension. Therefore, any weaker parameter 
than a ULAC function would have to fail. 

By the converse of the Hahn-Mazurkiewicz Theorem, when Z) is a bounded and 
simply connected domain, dD is locally connected if and only if it is the image of 
a continuous map on [0, 1]. However, it follows from the results in [4] that a ULAC 
function for dD can not be computed from such a map. Thus, a ULAC function 
appears to be an optimal parameter for the computation of boundary extensions of 
conformal maps. 
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The paper is organized as follows. In Section [3l we prove the complex analysis 
results which underpin our approximations to the values of boundary extensions. 
In Sectional we further discuss our approximation schemes including the approx- 
imations of arcs. Finally, in Section [5l we present our algorithm for computing 
boundary extensions. 

2. Preliminaries 

When S CAx B,\et 

dom(S') — {x : 3y £ S} 

ran(5) = {y : 3x {x,y) e S}. 

Let Dr(z) denote the open disk with center z and radius r. When X C C, let 

Dr{X) = U Dr{z). 

When X,Y C C are non-empty, let 

d{X,Y) ^M{\z-w\ : z £ X A w £Y}. 

Let d{p,X) — d{{p},X) when p e C and X C C is non-empty. When X C C is 
bounded and non-empty, let 

diam(X) = sup{|z — w| ; z,w € X}. 

When X C C and z e X, let Cz{X) denote the connected component of z in X. 

A chain is a finite sequence of sets (t/i, . . . , C/„) such that Ui n C/j 7^ whenever 
|i — j'l = 1. It is a simple chain if Ui ClUj $ precisely when |j — j| = 1. 

A continuum is a compact and connected metric space. A Peano continuum is 
a locally connected continuum. 

We will need the following results. Proofs can be found in Chapter 3 of [8]. 

Theorem 2.1. //{C/q}qg/ *s covering of a connected space X by open sets, and 
if ^iV G ^! then there exist ai,...,afc € / such that (Uai, ■ ■ ■ ,Ua^) is a simple 
chain from x to y. 

Theorem 2.2. Let X be a topological space. Then, X is locally connected if and 
only if for every open U C X, each connected component of U is open in X. 

Theorem 2.3. Every connected open subset of a Peano continuum is arcwise con- 
nected. 

3. Estimation results 

Let Cr = Dr{l)r\B. 

Again, suppose we have a conformal map of D onto a bounded domain D. We 
begin by considering the sort of situation illustrated in Figure [TJ Namely, let 
Asg = (/)[Cso], and suppose we have an arc C from a point p G As„ to a point 
q G dD such that C D dD = {q}. Our goal is to show that if the arc traversed 
by as t ranges from 1 — sq to 1, stays far enough away from the arc C long 
enough, then it will never cross C. We aim also to do this in a computable manner 
in that we exactly specify how far away this must be and for how long. 
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Figure 1. 



When 0< so,r <1 and Nq > 0, let 

27riVo 



x{so,No,r) = - 

In(so) - ln(l - r) 

m{soNo,r) = 2^/x{soNo,r) 

Note that m{so, No,r) > if sq > 1 — Note also that m{sQ, No,r) 0+ when 
r 1~. 

Theorem 3.1. Let (p he a conformal map of 1} onto a bounded domain D <Z C 
Let ro, So be numbers such that 0<1 — so<ro<l, and let Nq be larger than the 
area of D. Let Ag^ = (^[C^j,]. Suppose C is an arc from a, point p € As„ to a point 
q e dD such that C n dD = {q} and such that d{4>{t),C) > m{so, No,ro) for all 
te[l- so, ro]. Then, (l){t) ^ C for all i e [1 - so, !]• 

Proof. By way of contradiction, suppose otherwise. We use a Length-Area argu- 
ment. Let to be the smallest number such that (/>(to) € C. Thus, to > ro- Let 
a = (j)~^[C n D]. Hence, a is an arc in D from 4'~^{p) S Cgg to to S [?'o,l]- It 
follows that dDr{l) (la ^9 whenever 1 — ro < r < sq. 

Let 7(r, 6*) = 1 + re*^. Wc claim that if 1 — ro < r < Sq, and if ■y{r, 9) € cr, then 
^{r,0) ^ R. For, suppose otherwise. Since \j{r,9) — 1| = r, ^{r,0) = 1+r. Since 
a CI}, 7(r, 0) = 1 — r. But, since 1 — ro < r < sq, 1 — so < 1 — r < ro and so, by 
assumption, — r) ^ C. But, C C by definition of cr. 

So, for each r S [1 — ro,so], let 9r be the number such that Or G [7r/2, 37r/2], 
7(r, ^r) € a, and the distance from 7(r, Or) to the ar-axis is minimal among all such 
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points in a Ci Dr{l)- Hence, 6r tt- Without loss of generality, assume 6r > tt 
whenever re [1 — tq, sq]. 

Abbreviate to(so-^Oj''o) by m and x(soM),?'o) by x- For each r G [1 — vq^sq], 
7(r,7r) e M. So, 

(3.1) m < |(/.7(r,^)-07(r,0OI- 

At the same time, the right side of p.ip is not larger than the length of the arc 
traced by (/)7(r, 6) as 9 ranges from tt to 9r- The length of this arc is 



d_ 

89 



Mr, 9) 



d9 



\<j>'j{r, 9)\rde. 



de 



So, by the Schwarz Integral Inequality, 



\<j>'j{r,9)\r d9 



< 



|0'7(r,0)|2 de / d9 



= r\9r-7r) / |(/.'7(r,0)|2rf0 



< 27rr2 / \(j)'-f{r,9)\^d9. 



Therefore, 



/ -dr < 

27r Ji^ro r 



so /.e 



l-ro Jtt 



|0'7(r, 9)\^rdedr. 



Since r is the Jacobian of 7, 

(-so r^r 



b'j{r,9)\^rd9dr < 



dA. 



Di(i)nD 



However, by the Lusin Area Integral (see e.g., Lemma 13.1.2, page 386, of |6|), the 
latter integral is the area of (j)[Di(l) n D] which is smaller than Nq. So, it follows 
that 



f ° 1 2tt 
In(so) - ln(l -ro)= - dr < Nq^. 



Hence, 



27riVo 



In(so) - ln(l - ro) 



= X- 



But, by definition, > x- This is a contraction, and the proof is complete. □ 



When J is a Jordan curve, let Int(J) denote its interior. 

Corollary 3.2. Suppose 4>, D, sq, tq, and Ag^ are as in Theorem \ 3.1l Suppose 
Cj , pj , qj , and a are as in Figure That is, 
(1) (ji,(j2,0(l) are distinct, 
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(2) Cj is an arc from pj G As„ to qj e dD such that Cj n (t U dD) — {pj, qj} 
where t is the subarc of Ag^ from pi to p2, 

(3) a C dD is an arc from qi to q2 such that (piro) € Int(o' U Ci U C2 U r), 

(4) (j){so) is between pi and p2 on As^, and 

(5) d(0(i),Ci UC2) > m(so,^o,r-o) for all t e [l-so,ro]. 
Then, G Int(cr U Ci U C2 U r). 




dD 



Figure 2. 

Theorem 3.3. Suppose <j) is a conformal map of D onto a bounded and simply 
connected domain D. As in Theorem \3.1[ let As = 4>[Cs\ for all < s < 1. Then, 
for each e > 0, there exist SQ,rQ, Cj, pj, qj, a, and t that satisfy the hypotheses of 
Corollaru \S.2\ and such that the diameter 0/ cr U Ci U C2 U r is smaller than e. 

Proof. Let Sq, li, h, etc. be as in Figure [31 Let pj — 4>{p'j), and let qj — 4>[q'j). Let 
Cj — (t>[lj]. Let t' be the subarc of dCso from p'^ to p'2 that contains 1 — sq, and let 
T = 4)[t']. Thus © and (g]) hold. 

We claim that for each i5 > 0, there exist li, I2 as shown so that d{li,l2) < S 
and so that qi, q2, and 0(1) are distinct. For, suppose otherwise. It follows that 
there is a subarc a of i9D such that (j)\a\ = {^(1)}. Let ai C D be an arc that 
intersects 9D only at the endpoints of a. Then, a U ai and (/)[a U ai] are Jordan 
curves, and (j) maps the interior of the former onto the interior of the latter. Thus, 
by the Caratheodory Theorem, (p is injective on a U ai- a contradiction. So, we 
assume d{li,l2) < S and that ^ holds. 

Since D is simply connected, the boundary of D contains at most two arcs from 
qi to 92- Exactly one such arc has the property that the interior of the Jordan 
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P' — 



/ 1 

1^' q' 



Figure 3. 



curve obtained by adding it to Ci U C2 U r contains (j>{ro). Label this arc a. Let 
(3 be the subarc of dD from q[ to (72 that contains 1. It foUows that a C (j>[/3]. Let 
J = cr u Ci U C2 U T, and let Ji = r' U U ^2 U /3. Thus, J C (j)[Ji]. Since is 
uniformly continuous on D, the diameter of (/)[Ji] tends to as (sq, 5) — ^ (0, 0). So, 
we can assume ^ also holds and that the diameter of cr U Ci U C2 U r is smaller 
than e. 

Now, keeping sq, h, and h fixed, let tq — >■ 1~. Thus, m(so, A^o, ''0) 0"*". 
Whenever i g [1 — sq: 1]: 0(0 ^ Ci U C2. It follows that we can additionally choose 
a number tq so that (O holds. This choice does not affect any of ([Ij - ([4]). □ 



4. Approximation schemes 

Whenever p e C is approximated by a rational rectangle i?, the diameter of R 
shall be referred to as the maximum possible error of the approximation. 

Whenever a compact set X C C is approximated by a tight cover . . . , i?^}, 
the maximum possible error of this approximation is defined to be 
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max{diam(i?i), . . . , diani(i?fc)}. We also refer to this quantity as the diameter of 
the cover. 

Whenever we approximate a continuous function by a finite collection of pairs 
as in Section [T] we do not quantify the error of such an approximation. However, 
we say that one such approximation is no worse than another if every function 
approximated by the first approximation is also approximated by the second. 

In addition to points, compact sets, and continuous functions on D, we will 
need to approximate arcs. We will need to consider two situations: arcs through 
C and arcs through dD. When dD is bounded and locally connected, it is a 
Peano continuum. Thus, we also consider how to approximate arcs through Peano 
continua. 

Arcs are of course compact sets. However, to approximate them by tight covers 
does not yield sufficient information for most applications. For example, J. Miller 
has shown that there is an arc A C C that is computable as a compact set but 
which has no computable injective parameterization [12] . Hence, we turn to the 
following definitions which are essentially from [5]. 

A witnessing chain is a finite sequence (m, _Ri, . . . , Rk) where m g N and (-Ri, . . . , Rk) 
is chain of rational rectangles. 

However, when X is a planar Peano continuum and 5 is a ULAC function for 
X, then a witnessing chain for X restricted by g is a finite sequence of the form 
(m, . . . , Rk) such that each Rj is a rational rectangle, Ci X, . . . , Rk Ci X) is 
a chain, and diam(i?j) < 2"^*^™^ whenever 1 < j < k. (Thus a witnessing chain for 
X restricted by 5 is a witnessing chain.) 

We now make some notation. Let uj = (to, _Ri, . . . , Rk) be a witnessing chain. 
We let: 

k 

14 = [jD,-^{R,) 

i=i 
rriui = m 

k^ — k 
Ru.j = Rj 

An arc chain is a finite sequence of witnessing chains (wi,...,w„) such that 
(Kji , • ■ • , VL„) is a simple chain. 

But, when X is a planar Peano continuum for which we have a ULAC function 
5, then an arc chain for X restricted by g is a, finite sequence of witnessing chains 
for X restricted by g, (wi, . . . , a;„), such that (K;^ , . . . , VL„) is a simple chain and 
such that X n R^j^k^. n -R(^j_|_i,i 7^ whenever 1 < j < n. 

We make some notation. Let p — (wi, ... be an arc chain. We let: 

- I 

yp.3 = 

The diameter of p is then defined to be maxj diam(Vp j). We denote this by diam(p). 
The sets Vp j are called the links of p. 
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Suppose A is an arc, and that {Ui, . . . , C/„) is a simple chain. We say that A goes 
straight through {Ui, . . . , Un) if one endpoint of A belongs to Ui, the other belongs 
to Un, and there do not exist pi,p2,P3 G A and I < j < n such that p2 is between 
pi and p3 on A, pi^ps € Uj — Uj+i and p2 € — Uj. i.e. A does not "backtrack" 
as it goes through (C/i, . . . , [/„). 

We approximate an arc A C C by an arc chain (wi, . . . ,a;„) such that A goes 
straight through (VLi , ■ ■ ■ , VLj„)- 

However, whenever A is contained in a planar Peano continuum X for which we 
have a ULAC function g, we may approximate A by an arc chain for X that is 
restricted by g. In either case, the maximum possible error of the approximation is 
defined to be the diameter of the arc chain. 

The results in the remainder of this section are adaptations of results in [5] . 

Lemma 4.1. Let X be a planar Peano continuum. Suppose u ~ (m, . . . ,-Rfc) 
is a witnessing chain for X restricted by g, and that 1 < j < k. Then, for all 
z,w e Rj n X, 

CM-^iRj) nx)^ CM-^4Rj) n x). 

Proof. Let U = D^-^iyRj). Since diam(i?j) < 2-f(™), it follows that w & Xn 
D2-g(m){z). Since D2-m{z) C U, it follows that 

c,{D2-r^{z)nx) c c^unx). 

Since g is a ULAC function for X, Xn i^a-sc™) i^) ^ Cz{D2-m (z) n X). It then fol- 
lows that w G Cz{Ur\X). So, Cz{Ui^X) and Cw{U CiX) are connected components 
of U n X having a point in common, namely w. Hence, they are equal. □ 

We make some more notation. Suppose a; is a witnessing chain for a planar 
Peano continuum X and is restricted by g. We let 

for any z G R^j n X. By Lemma [4.11 the choice of z is inconsequential. We then 
let 

k 

Let p — (wi, . . . ,uji) be an arc chain. We let: 

3 

Theorem 4.2. (1) If p is an arc chain, then there is an arc approximated by 

(2) // p is an arc chain for X restricted by g, then there an arc A Q X that 
is approximated by p. Furthermore, if pi ^ X belongs to a rectangle in the 
first witnessing chain in p, and if P2 ^ X belongs to a rectangle in the last 
witnessing chain in p, then there is an arc A <Z X from pi to p2 that is 
approximated by p . 

Proof. We first prove ([2]). The proof of ([T]) is a straightforward simplification. 

Let p = (wi, . . . ,uJi). Thus, pi S Cui^ and p2 S Coji. It follows from Theorem [221 
that each Cui^ is a connected and relatively open subset of X. Hence, by Theorem 



10 



TIMOTHY H. MCNICHOLL 



12.31 they are all arcwise connected. It also follows from the definition of 'arc chain 
for X restricted by g' that {C^ji , • • ■ , Cuj, ) is a simple chain. Let qj e C^j^ H C^^^^j 
whenever 1 < j < k. 

C^-^ contains an arc Ai from pi to qi. Whenever \ < j < k ~ \, C^-^-^ contains 
an arc from qj to qj+i- C'l^^. contains an arc Ak from qk-i to p2. 

Let q'l be the point in Ai D A2 such that there is no point of A2 between pi and 
q'l on Ai. Let Bi be the subarc of Ai from pi to q[. Whenever 1 < j < k, let q'j 
be the point in Aj D Aj+i such that there is no point of Aj+i between q'j^i and 
on Aj. Let Bj be the subarc of Aj from qj_-^ to q'j. Therefore, H Bj = {q'j^i} 
and n {Aj+i U . . . U Afc) = {g^}. 

Let Bf^ be the subarc of A^ from g^ j^ to p2- We may then take A = [jj Bj. □ 

Suppose u! — (to, R, Ri, . . . , Rk) is a witnessing chain. If p e and if g £ Rt, 
then we say that a; is a witnessing chain from p to q. 

Lemma 4.3. Suppose g is a ULAC function for a planar Peano continuum X . Let 
U ^ C be an open set, and suppose z,w are distinct points in the same connected 
component of U H X . Then, there is a witnessing chain uj from z to w such that 
QU and such that uj is a witnessing chain for X restricted by g. 

Proof. It follows from Theorems l2.2l and l2.3l that there is an arc A from ztow which 
is contained in U D X. There exists m such that D2-m{A) C U. For each p E A, 
there is a rational box Sp such that p G Sp, £)2-m (Sp) C U, and diam(5'p) < 2~^(™\ 
By Theorem 12.11 there exist pi, . . . ,pt £ A such that {A n Sp^ , . . . , A Ci SpJ is a 
simple chain. Set uj — (m, Sp-^, . . . , Sp^). It follows that cj is a witnessing chain 
from z to w. It also follows that K, C [/ and that a; is a witnessing chain for X 
restricted by g. □ 

Theorem 4.4. Suppose A C_ C is an arc and e > 0. 

(1) There is an arc chain p that approximates A with error at most e. 

(2) // A is contained in a planar Peano continuum X , and if g is a ULAC 
function for X , then there is an arc chain p for X restricted by g such that 
p approximates A with error at most e. 

Proof. We prove the second claim first. The proof of the first claim is a straight- 
forward simplification. 

Divide A into n subarcs Ai, . . . , A„ of diameter at most e/3. Label these so that 
Aj^nAj^ ^ precisely when \ji—j2\ — 1- Let < S < min{d{Aj^, Aj^) : |ji— ^2! > 
1}. Cover Aj by finitely many open disks of radius 61 where Si — min{e/3, ^/2}; 
let Uj be the union of these disks. Then, (C7i, . . . , C/„) is a simple chain and A goes 
straight through it. It follows from Lemma [4.31 that there exist witnessing chains 
ioi X u!i, . . . ,uJn that are restricted by g and such that Aj C \{j ■ C Uj . The proof 
is complete. □ 

A circular chain is a finite sequence of witnessing chains (wi, . . . ,u!k) such that 
Vuj^^ n VLj2 if and only if ji - j2 = +1 mod k. Again, the sets VLi , ■ ■ ■ , Ki^ 
are called its links. For all j G N, let j denote the smallest positive integer that is 
equivalent to j modulo k. We say that an arc A goes straight through a circular 
chain (wi, . . . ,LUk) if there exist ji < . . .j„ such that {^jJj^, ■ ■ ■ j'*^) is an arc chain 
that A goes straight through. We then say that a Jordan curve goes straight 
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through (wi, . . . jWfc) if aU of its subarcs go straight through (cji, . . . ,LOk)- Again, 
the diameter of a circular chain is defined to be the largest diameter of its links. 

We approximate a Jordan curve J by a circular chain (wi, . . . ,a;„) such that J 
goes straight through (Kj^, . . . , V^^). The maximum possible error is defined to be 
the diameter of the circular chain. By elaborating slightly on the proofs of Theorems 
I4.2l and l4.4l it can be shown that every Jordan curve is approximated with arbitrary 
precision by circular arc chains and that each circular chain approximates at least 
one Jordan curve. 

In any case, if we map say '^^J into V^^- , then we obtain an approximation 
to a parameterization of A/ J. 

We use the notation X to denote an approximation to X. 

An approximation to an object typically approximates many objects. For exam- 
ple, a rational rectangle approximates each point it contains. Accordingly, we let 
{X) denote the set of all objects approximated by X. 



5. Algorithms 
5.1. Algorithm 1. We are given the following as input: 

• An approximation </> to a conformal map </> of D onto a bounded domain D 
whose boundary is locally connected. 

• An approximation dD to the boundary of D. 

• An approximation 5 to a ULAC function for dD, g. 

It is required to produce an approximation to 0(1). 

The strategy of the algorithm is to produce approximations to arcs Ci , C2 1 
(T so that there exist Cj € (Cj), t G (t), a G ((?) and rational numbers so,rQ 
such that the hypotheses of Corollary [32] hold. To this end, we introduce a notion 
of substantiation as follows. That is, what it means for the input data {dD,(j3,g) 
to substantiate an approximation to some object. Intuitively, this means that the 
input data contain sufficient information to determine that the intended approxi- 
mation to the given object really is an approximation to that object. 

When we write that {dD,(j),g) substantiates (j){z) we mean that there is a pair 
{U, V)€4) such that z S C/ and ^) = V. 

When we write that {dD,(j},g) substantiates (j)[s,r], we mean that there are 
pairs {Ui, Vi), . . . , (?7„, Vn) £ <j) such that {Ui, . . . , C/„} tightly covers [s, r], and 

^]^{Vu...,Vn}. _ _ 

When we write that {dD, (j),g) substantiates Ci, we mean that Ci is an arc 
chain (wi, . . . , Wfc), each rectangle in ujk belongs to ran((^) U dD, and that for each 
1 < n < k, each rational rectangle in LUn belongs to ran((/)). 

When we write that {dD,(f),'g) substantiates C2, we mean that C2 is an arc 
chain (wi, . . . , cjfe), each rectangle in cji belongs to ran((/)) U dD, and that for each 
1 < n < k, each rational rectangle in ujn belongs to ran((/)). 

We say that r € Q is bounded by m £ N if there are n,d £ {—rn, . . . , m} such that 
r = n/d. We say that a rational point p is bounded by m if each of its coordinates 
is bounded by m. We say that a rational rectangle R is bounded by m if each of 
its vertexes is bounded by m. 
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Suppose oj = {m, . . . , Rk)- We say that ((/>, dD, g) substantiates that ui is a 
witnessing chain for dD restricted by g if for every 1 < j < fc i?j fl Rj+i contains a 
rectangle in dD and m G dom(^) and diam(ii'j) < 2^^^^™) for every 1 < j < k. 

Let us say that {(j),dD,g) substantiates a if it substantiates that each wgj is a 
witnessing chain for dD restricted by g and for each 1 < j < h 

contains a rectangle in dD. 

Let fco = ri/diam(aD)]. Let so e (0, 1) n Q. Fbc 1 < fc < fco- Divide Cs^ into 

2k arcs of equal length. Label these in clockwise order as X-k, • • • , A_i, Ai, . . . Afe so 
that those with negative indices lie below the x-axis and those with positive index 
lie above the a;-axis. When we write that {(f), dD, sq) substantiates r, we mean that 
we can choose k < ko and 1 < fc' < fc so that for each j = —k',...,k' we can 
additionally choose (C/p\ V"/-'-'), . . . , {Un,\Vn{'') G (f) such that A^ is tightly covered 
by : s = 1, . . . , rij} and the rectangles in the j-th witnessing chain of r are 

{vl^\...Mi^}- 

Our algorithm proceeds as follows. First, set 

A^o = 7rdiam(|JaD)2. 

Thus, A^o is larger than the area of D. For each 1 < ^2 < /ci < fco, we do 
the following. First, set sq = and ro = 1 — Wc then search for d, Ci, C2, 

(t>{l — So), and (t>[l — so, ro] that are substantiated by {dD, (j),g) and r substantiated 
by {dD, 4>,g, sq) so that the following hold. 

(1) Each witnessing chain in a, Cj, f has the property that each of its rational 
rectangles is bounded by fco- (This condition ensures that the search is 
bounded.) 

(2) The last witnessing chain of C\ contains one of the rectangles in the first 
witnessing chain of a. 

(3) The first witnessing chain of C2 contains one of the rectangles in the last 
witnessing chain of a. 

(4) The first witnessing chain of C\ contains one of the rectangles in the last 
witnessing chain of t. 

(5) The last witnessing chain of C2 contains one of the rectangles in the first 
witnessing chain of f . 

(6) — So) is contained in an intermediate link of r. 

(7) J =dj (Ci,(T, 6*2,?) is a circular chain. 

(8) There is an integer < t < fco and U\,...,Un G dom(0) such that: 

• diam((T) < 2~* 

• n CT, (72 n (t) < 2-f(*) 

• Ui, . . . ,Un & dom((/)) cover [0, 1 — 2so] 

• Uj n [1 - So, 1] = 

2"* + diam(Ci U (72 U f) < diam(|J Uj). 

j 

(These conditions together imply that (?i(ro) C Int( J) for any J G ( J).) 
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(9) d{[jcp[l- so,ro],{jCiU[jC2) > m{so,No,ro). 

If none are found, then output (— [diani(9Z3)] , [diam(9D)] )^. Otherwise, let R^^ g 
label the rectangle 

(a - 2"™, b + 2-™) X (c - 2-™, d + 2-™) 

where (f>{ro) ~ (a, 6) x (c, d) and m is the largest number such that diani(lj J) < 
2-™. We then output 

(5.1) n R 

sa,ra,C\,Ci.a .T 

We now demonstrate the correctness of the algorithm. That is, we show that if 
R is output, then 0(1) € R. This is clear if no sqj ^Oi e^c. are found. So, suppose 
(|5.ip is output. It suffices to show that 0(1) belongs to each R^^^ raChCo.'S ?• 
this end, it suffices to show that there are arcs a e (ct), r G (r), Ci € (C'l), and 
C2 S (C2) as well as points Vj^qj such that the hypotheses of Corollary 13.21 hold. 

By construction, (r) contains a subarc of A^^. Let t' denote the largest such 
subarc. 

The first witnessing chain of r contains a rational rectangle that also belongs to 
the last witnessing chain of C2. This rectangle contains a point of t'. It similarly 
follows that there is a rectangle in the first witnessing chain of C\ that contains a 
point of t'. It now follows that there is an arc C( G (C'l) from a point ])\ on r' to 
a point on (9£). C( can be chosen so that belongs to the last link of r and 
belongs to the first link of ct. It also follows that there is an arc C'^ S (C2) from a 
point ^2 S 91) to a point p'2 ^ t' . Furthermore, C2 can be chosen so that 93 belongs 
to the last link of a and p'2 belongs to the last link of r. Let pi be the last point 
on C[ that belongs to r'. Hence, pi belongs to a rectangle in the last witnessing 
chain of r. Let p2 be the first point on that belongs to r'. Hence, p2 belongs to 
a rectangle in the first witnessing chain of t. 

There is an arc a' € ia) from q'^ to Hence, cr' C dD. Let gi be the first point 
on Ci that belongs to cr', and let 92 be the last point on C2 that belongs to a' . 
Thus, qi belongs to the first link of tr, and 92 belongs to the last link of a. 

Let cr be the subarc of cr' from qi to (72- Let r be the subarc of r' from pi to p2- 
Let Ci be the subarc of C( from pi to qi . Let C2 be the subarc of from p2 to 
(72- It follows that CT, T, Cj, pj, satisfy the hypotheses of Corollarv l3.2l 

We now demonstrate the convergence of the algorithm. That is, we show that if 
R is any rational rectangle that contains 0(1), then there exist 0, dD^ g such that 
when the algorithm is provided (0, dD, g) as input, it produces a rational rectangle 
R' that is contained in R. To this end, choose ei > such that Z3ej(0(l)) C R. By 
Theorem 13.31 we can choose pj, qj, rg, sq, a, t, Ci, and C2 so that the hypotheses 
of Corollary 13.21 hold and so that the diameter of cr U r U Ci U C2 is smaller than 
ei/3. _ _ 

Let g be a ULAC function for dD. We first choose a, r, Ci, C2, and g as follows. 
To begin, we choose these arc chains so that they have diameter smaller than ei/3. 
We can also choose them so that J =df (t, Ci , ct, C2) is a circular chain. Let g be the 
smallest initial segment of g used in the formation of ct. It is now possible to choose 
and dD so that (0, dD,g) substantiates Cj, a, and (0, dD,g, sq) substantiates f. 
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Note that the diameter of IJ J is smaUer than ei. We can choose 4> so that 

(4>,dD,g) substantiates an esitmate ^(ro) whose diameter is smaller than ei. It 
now follows that if ei is sufficiently small, then R ~ ^ ^ ~ d R. Convergence 
is thus demonstrated and the verification is complete. We note that we can also 
choose 0, dD, g so that they improve on any other approximations to </>, dD, g. 

We also note that the algorithm has a monotonicity property^ That is if (/>, 
9D, g are approximations to dD^ g respectively, and if (pi, dDi, gi are also 
approximations to (j), dD, g, then the output when {<p, dD, g) is provided as input 
is no worse an approximation to 0(1) than when {(j)i,dD\,g{) is provided as input 
provided that the first set of approximations is no worse than the first. 

5.2. Algorithm 2. Suppose we are given 0, dD, g, and p as input and that 

• dD approximates the boundary of a bounded, simply connected domain D, 

• g approximates a ULAC function for g, 

• 4> approximates a conformal map of P onto D, and 

• p approximates a point p G dH. 

It is required to produce an approximation to 0(p). 

If the diameter of pTlcJID) is not smaller than v^, then output (— [diam(3£')] , [diam(9D)] )^. 
Otherwise, compute ai , a2 such that 

pn dB = {e'^ : ai < 6 < 

and so that a2 — ai < 7r/2. Let 

So, if p = e'^ and if tp{z) ~ 0(e*^z), then ip e (■!/>). We then run Algorithm 1 on 
{i^,db,g). 

5.3. Algorithm 3. Suppose we are given 0, dD, g, p such that 

• dD approximates the boundary of a bounded, simply connected domain D, 

• 'g approximates a ULAC function for dD, g, 

• (j) approximates a conformal map of D onto D, </>, and 

• p approximates a point p e D. 

It is required to produce an approximation to <p{p), <p{p)- We will take advantage 
of the following observation: the rectangle produced by Algorithm 1 contains not 
merely but also 0[ro. 1]. 

We first set Ro = (-[diam(a5)l, rdiam(aD)])2. Hence, D C Rq. 

We then consider the case where p contains no point of d'D. If the closure of p 
contains a point of 9D, or if p is not contained in any rational Carleson rectangle 
in dom(0), then we output Ro- Otherwise, we output 



f]{V I 3U iU,V) e A U}. 



Suppose on the other hand that p (1 dV) ^ 0. We then run Algorithm 2. But, 
when running Algorithm 1 as a result, we restrict ro to values smaller than 

min{Re(e' ^z) : z ep A -a2 <9 < -ai}. 

Correctness and convergence follow immediately. 
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